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Assessing  friction  stress  on  a  liquid  lubricant  by  stable  isotope  analysis 
1.  Introduction 

Wear  and  tear  of  moving  meehanieal  parts  is  eostly  to  individuals,  eompanies  and 
nations,  sealing  from  the  thousands  to  billions  of  dollars  in  annual  loss  [1].  Lubrieants  are 
the  front  line  of  defense  against  meehanieal  failure.  Therefore,  most  large  organizations 
develop  lubrieation  management  plans  to  ensure  equipment  reliability,  reduee  down-time, 
and  minimize  both  waste  generation  and  energy  use  [2],  A  proaetive  strategy,  referred  to 
as  eondition-based  maintenanee,  involves  performing  preventative  maintenanee  only 
when  operating  eonditions  warrant  (e.g.  [3]).  Implieit  in  this  strategy  is  the  availability  of 
a  range  of  independent  methods  for  analyzing  and  monitoring  maehine  and  lubrieant 
eondition. 

Lubrieant  degradation  analysis  generally  follows  two  basie  strategies;  (1)  determine 
the  distribution  pattern  of  eomponents  of  the  lubrieant  itself  (lubrieant  and  additives),  or 
(2)  determine  the  eontaminants  from  wear  within  the  lubrieant  matrix  (e.g.  metals,  fuel, 
moisture,  ete.).  The  latter  may  provide  a  hardware  eondition  snapshot  and  is  by  far  the 
most  eommon  test  in  terms  of  both  historieal  and  present  use.  In  addition,  most  sensor 
development  teehnologies  revolve  around  assaying  wear  eomponents  and  eontaminants 
in  the  lubrieating  oil  matrix  {c.f.  [4]). 

Lubrieant  eomponent  analyses  may  be  less  intuitive  to  interpret  and  assess.  However, 
it  may  be  most  applieable  to  determining  the  lubrieant' s  or  additive's  useful  life  and  may 
provide  information  on  degradation  rates  and  states.  Teehniques  have  been  developed  to 
assess  struetural  ehanges  to  lubrieant  and  additive  ehemistries  or  eomponent  evolution. 
The  most  widely  applied  teehniques  have  been  infrared  (IR)  absorption,  Raman  seattering 
and  ultraviolet  fluoreseenee  [5].  Advanees  in  both  IR  [5-8]  and  fluoreseenee  speetroseopy 
[5]  eoupled  with  multivariate  statisties  have  enabled  lubrieant  breakdown  traeking  from 
modeled  frietion  stress.  Tribologieal  investigations  using  speetral  teehniques  have 
allowed  eorrelations  between  frietion  stress  (engine  use,  heat,  ete.)  and  eomponent 
oxidation,  deeomposition  of  ester  funetionality,  additive  breakdown,  and  breakdown 
produet  emergenee  [5,  8-10].  While  eaeh  of  the  aforementioned  methods  has  strengths, 
additional  indiees  whieh  ean  relate  wear  (under  varying  eonditions)  to  lubrieant  speetral 
eharaeteristies  would  be  of  both  eeonomie  and  seientifie  importanee. 

In  meehanieal  systems,  energy  may  be  eonverted  into  frietion  where  parts  move  in 
relation  to  one  another.  This  frietional  energy  ean  then  be  transferred  to  a  lubrieant  under 
all  fluid  fdm  lubrieation  regimes  {i.e.  hydrostatie,  hydrodynamie,  boundary,  ete.)  [11]. 
Frietion  and  assoeiated  heat  may  be  suffieient  to  initiate  lubrieant  deeomposition.  If 
deeomposition  involves  the  breaking  of  bonds  (e.g.  oxidation,  ester  eleavage,  ete.), 
reaetions  have  the  potential  to  involve  kinetie  isotope  effeets  (KIE).  For  ehemieal 
reaetions,  the  mass  differenee  of  heavier  isotopes  relative  to  lighter  isotopes  of  the  same 
element  neeessitates  higher  initial  aetivation  energies  (for  review  see  [12]).  For  instanee, 
there  is  an  8%  mass  differenee  between  the  heavy  earbon  isotope  (  C)  relative  to  the 
light  isotope  ('^C).  The  heavy  isotope  of  hydrogen  (deuterium,  or  ^H)  is  100%  heavier 
than  the  light  isotope:  ^H.  In  a  ehemieal  reaetion,  this  phenomenon  may  lead  to 
isotopieally-lighter  produets  and  an  isotopieally-heavier  residual  pool  of  initial  eompound 
during  degradation  [12]. 
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While  originally  employed  in  enzymology,  kinetie  isotope  effects  have  recently  been 
used  in  “field”  studies,  principally  to  assess  biological  contaminant  degradation  [13].  In 
these  studies,  if  the  “original”  contaminant  pool  becomes  isotopically-heavier  as  its 
concentration  decreases,  biological  degradation  {i.e.  bond  cleavage)  is  inferred. 
Additionally,  using  multivariate  chemometric  methods  to  assess  isotopic  variation  in 
multi-component  compound  mixtures,  one  can  determine  if  two  sources  are  or  are  not 
statistically  related  [14].  If  the  transfer  of  frictional  energy  to  lubricants  and  subsequent 
oxidation  or  bond-cleavage  imparts  an  isotope  effect,  isotopic  shifts  in  lubricant 
components  could  serve  as  a  friction  exposure  diagnostic.  To  our  knowledge,  this 
hypothesis  has  not  been  tested.  We  thus  hypothesize  that  lubricant  component  and 
degradation  product  hydrogen  and  carbon  isotope  ratios  will  be  impacted  by  friction 
stress  and  that  multivariate  statistical  models  can  be  applied  to  allow  an  independent 
index  of  lubricant  sliding  exposure,  perhaps  providing  a  novel  test  mechanism  and 
allowing  independent  validation  of  other  lubricant  wear  tests. 

2.  Methods 

2.1.  Lubricant  tests 

A  pentaerythritol  ester  lubricant  (Hercolube-A,  Hercules,  Inc)  was  chosen  which  has 
a  widespread  military  use  (aviation)  and  demonstrated  readily  identifiable  peaks  via 
standard  temperature  gas  chromatography  (GC)  analysis.  Hercolube-A  (Here -A)  was 
diluted  at  a  1:100  (v/v)  ratio  with  dichloromethane  for  GC  and  FTIR  analysis.  A  pin-on- 
disk  tribometer  (CSEM,  Neuchatel,  Switzerland)  was  used  for  exposing  the  lubricant  to 
friction  stress.  A  10  pL  neat  aliquot  of  Here -A  was  pipetted  through  the  center  of  the  ball 
holder  onto  the  slowly  rotating  disk  at  the  beginning  of  each  experiment  in  a  manner 
similar  to  Dugger  et  al  [15].  Both  ball  (1/4”  diameter)  and  disk  were  AISI  52100  steel, 
and  a  5  N  load  was  applied  for  the  duration  of  each  experiment.  Runs  were  made  in 
ambient  air  at  25±3  °C;  relative  humidity  was  kept  below  10%  by  addition  of  desiccant  in 
the  tribometer  chamber.  The  sliding  radius  was  fixed  at  15  mm  for  each  run,  and  at  100 
cycles/min  resulted  in  a  linear  sliding  speed  of  0.15  m  s'\  Triplicate  runs  were  conducted 
for  the  annotated  cycles  (Table  1).  After  completion  of  each  run  and  a  0  cycle  control,  the 
ball  and  disk  were  immersed  in  dichloromethane  and  sonicated  for  30  minutes  for 
thorough  lubricant  to  solvent  transfer.  The  Here -A  containing  solvent  was  then 
transferred  to  a  calibrated  tube  and  dried  down  to  1  mL  using  a  Zymark®  Turbo Vap 
system  supplied  with  N2  gas. 

2.1.  Carbon  GC-Quadrapole  MS  and  GC-Combustion-lsotope  Ratio  MS  analysis 

Here -A  was  analyzed  using  a  custom-configured  GC-Combustion-Isotope  Ratio  MS 
(GC-C-IRMS)  system  similar  to  one  previously  reported  [14].  The  configuration 
consisted  of  a  Hewlett  Packard  6890  GC  -  5973  quadrapole  MS  outfitted  with  a  250  pm 
ID  30  m  Supelco  SPB-05  capillary  column.  A  post  column  4-way  splitter  (Valeo) 
allowed  a  20:80  split  between  the  quadrapole  MS  and  the  IRMS,  respectively.  IRMS 
flow  was  routed  through  a  Finnigan  GC  combustion  interface,  which  in  turn  interfaced  to 
a  Finnigan  MAT  Delta  S  IRMS.  Injections  were  splitless  to  minimize  potential  isotope 
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fractionation  in  the  inlet  [16].  The  inlet  was  set  isothermally  at  280  °C  and  operated  at  35 
psi  pulsed  pressure  for  the  first  30  seeonds  of  eaeh  run.  Flow  was  set  for  2  mL  min’' 
throughout  eaeh  run.  Oven  operating  eonditions  were  as  follows:  150  °C  for  5.00  min, 
ramp  to  280  °C  at  16  °C  min'\  hold  10  min,  ramp  to  310  °C  at  5  °C  min,  hold  12  min. 
Eaeh  sample  was  spiked  with  an  internal  isotopie  standard  (n-tetraeosane,  -31.548  ±  0.23 
%o),  whieh  eluted  immediately  before  lubrieant  peaks  and  ealibrated  by  thermal 
eonversion  elemental  analyzer  (TCEA)-IRMS  against  NBS -22.  An  isotopie  referenee 
CO2  gas  (-22.805  ±  0.12  %o)  was  injected  in  triplicate  direetly  into  the  IRMS  source  at  the 
beginning  and  end  of  eaeh  run. 

2.2.  Hydrogen  GC  -Isotope  Ratio  MS  analysis 

Here -A  was  analyzed  in  as  close  as  possible  GC  eonditions  as  for  5  C  analysis  {see 
above).  Setup  for  eontinuous  flow  5D  measurements  was  as  previously  deseribed  [17]. 
Briefly,  combined  replieates  were  run  in  triplieate  on  a  Varian  3400  GC  eoupled  to  a 
reduetion  fumaee  held  at  1450  °C.  In  this  configuration,  the  pyrolysis  interfaee  converted 
hydroearbons  to  H2/HD  gas  whieh  passed  through  a  Nafion®  dryer  before  flowing  into  an 
open  split  interfaee  to  a  Einnigan  DelatPlus  XE  IRMS.  H2  referenee  gas  was  injeeted 
through  an  open  split  at  the  beginning  and  end  of  eaeh  run  for  ealibration.  Duplicate  runs 
were  made. 

2.2.  FTIR  analysis 

To  determine  if  ETIR  analysis  was  eomparable  with  CSIA  to  assess  lubrieant  friction 
stress,  eaeh  lubrieant  sample  was  analyzed  both  with  sliding  exposure  alone  and  with  the 
addition  of  two  eommon  fuels  (whieh  might  likely  eontaminate  lubrieants  under  real- 
world  operational  conditions).  Eor  the  fuel  eontamination  tests,  0.4%  (v/v)  diesel  or  JP8 
fuel  was  added  to  the  original  and  friction  stressed  Herculube  A.  The  "original"  wear 
series  (Table  1)  and  fuel-added  series  were  pipetted  into  a  Cap2  window  ETIR  eell  with  a 
0.1mm  spaeer  and  analyzed  using  a  Mattson  IR-101 10  ETIR  speetrometer.  Raw  data 
were  correeted  for  the  methylene  ehloride  solvent  (run  alone)  and  exported  to  Matlab® 
for  modeling  and  analysis. 

2.4.  Data  Analysis 

Isotope  ratio  values  were  ealeulated  and  output  in  spreadsheet  format  by  IsoDAT  NT 
software  (Thermo-Eisher).  Isotope  ratio  data  for  eaeh  identified  peak  were  eollated, 
transferred  to  Matlab®  software,  and  corrected  for  the  internal  isotopie  standard  (IIS)  n- 
tetraeosane  using  equation  I  for  carbon  and  equation  2  for  hydrogen: 


6%eak  + (-35.458 -S^'Ciis)  (I) 

5Dpeak  + (-22.792 -SDiis)  (2) 

Two  general  statistieal  tests  were  applied  to  triplieate  isotope  data.  Eirst,  in  order  to 
determine  if  there  were  signifieant  differences  between  eaeh  treatment,  a  stepwise 
eoupled  prineipal  components  analysis  -  multiple  analysis  of  variance  (PCA-MANOVA) 


3 


was  conducted  [14],  Second,  to  determine  if  wear  cyeles  eould  be  related  to  isotopie  data, 
a  PCA  model  was  performed  on  the  raw  eolumnar  data.  In  the  PCA  analysis,  data  were 
autosealed  by  dividing  eaeh  value  by  its  eolumn  standard  deviation.  Finally,  peak(s) 
whieh  appeared  in  ehromatograms  after  wear  exposure  were  assessed  for  isotopie  shifts. 

FTIR  data  were  analyzed  by  PCA  and  the  resultant  seores  were  related  to  wear 
eyeles.  The  loadings  from  the  PCA  models  were  also  plotted  against  wavenumber  to 
assign  struetural  information  to  regions  whieh  varied  in  response  to  wear  treatment. 

3.  Results  and  Discussion 

3.1.  Hydrogen  isotope  ratio  shifts 

Ten  major  peaks  were  identified  and  used  to  perform  statistieal  analyses.  Corrected 
5D  values  ranged  from  -163  to  -47  %o.  Statistieal  analysis  revealed  only  one  instanee  in 
whieh  there  was  a  signifieant  differenee  between  samples  (Table.  2).  This  differenee  was 
between  the  initial  starting  material  and  the  lowest  sliding  cycle  exposure  sample  (24,000 
eyeles).  This  appears  to  be  anomalous,  as  one  would  expeet  higher  levels  of  sliding 
exposure  to  also  show  signifieant  differenees.  The  presenee  of  earbon  isotopie  shifts 
during  wear  tests  (see  below)  but  not  hydrogen  isotope  shifts  indieated  that  C-H  bonds 
were  not  likely  involved  in  initial  lubricant  degradation.  This  eould  oecur  if  the  frietion- 
indueed  ester  degradation  behaved  similarly  to  hydroearbon  thermal  eraeking,  in  whieh 
C-C  bond  eleavage  oeeurs  preferentially  to  C-H  bond  eleavage.  Beeause  C-C  bond 
energies  are  '-'200-300  kJ  mof^  while  C-H  bond  energies  are  -'300-500  kJ  mof',  the 
weaker  bonds  should  break  preferentially.  Considering  molecular  orbital  theory,  the 
earbon  radieal  in  a  C-C  bond  break  would  have  higher  stability  than  a  hydrogen  radieal  in 
a  C-H  bond  break  [18]. 

3.2.  Carbon  isotope  ratio  shifts 

Twenty  seven  disereet  peaks  were  identified  in  the  isotope  ratio  ehromatogram 
eollection.  Of  these,  three  were  not  found  in  every  chromatogram,  thus  were  not  included 
in  the  PCA-MANOVA  analyses.  The  PCA-MANOVA  model  showed  signifieant 
differenees  between  eaeh  treatment  (0,  24,  48,  96,  192,  and  384k  wear  eyeles)  with 
probabilities  usually  less  than  0.01  (Table  3). 

In  order  to  observe  trends  and  deeonvolve  tabular  data,  a  PCA  model  was  ereated. 

The  initial  PCA  model  showed  inflection  at  3  and  6  components.  The  6  eomponent  model 
had  eigenvalues  greater  than  1 .0,  and  eaptured  85%  of  the  variability  in  the  dataset.  A 
two  eomponent  model  was  "suggested"  by  the  software.  No  diseemable  trends  were  seen 
in  any  PCA  seores  aside  from  PCI  (with  models  from  2  to  6  eomponents).  Average 
seores  on  PCI  (from  triplicate  analyses),  which  accounted  for  44%  of  the  variation  in  the 
eompiled  data,  displayed  an  exponential  inerease  over  sliding  eyeles  (Fig.  1).  The  fit  was 
signifieant  with  a  0.96  goodness  of  fit  statistie  (r^).  A  biplot  showing  seores  and  loadings 
demonstrated  that  certain  peaks  were  instrumental  in  driving  the  model  PCI  scores 
between  0  and  384k  eyeles  (Fig.  2).  Here -A  eomponents  at  retention  times  (indieated  by 
blue  squares)  in  region  "A"  and  region  "B"  (outlined  in  magenta  on  the  figure)  were 
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responsible  for  inereases  or  deereases  in  isotopie  ratios  between  native  and  frietion 
stressed  lubrieant. 

Sereening  Here -A  eomponents  at  wavelengths  highly  loaded  along  PCI  (both  positive 
and  negative)  showed  eomponents  (at  retention  times  shown  in  figures)  whose  isotopie 
ratios  inereased  (Fig.  3)  and  deereased  (Fig.  4)  as  a  funetion  of  sliding  eyeles.  In  order  to 
be  eonsistent  with  the  kinetie  isotope  model,  an  inerease  in  any  residual  lubrieant  peak 
isotope  ratio  would  be  expeeted  {i.e.  isotope  ratios  of  peaks  that  deerease  in  abundanee 
due  to  frietion-indueed  degradation  should  beeome  more  positive  in  value).  Under  the 
relatively  mild  frietion  eonditions  applied  to  Hero -A,  one  “new”  peak  was  observed  after 
24K  eyeles.  The  5  C  values  for  this  peak  deereased  logarithmioally  in  relation  to  wear 
eyeles  (see  peak  944  in  Fig.  4;  =  0.96,  P  <  0.01).  A  linear  relation  also  had  a  relatively 

high  goodness  of  fit  (r^  =  0.91,  P  <  0.01). 

To  assess  overall  oonoentration  ohanges  during  frietion  stress,  peak  areas  were 
normalized  to  the  surrogate  standard  (tetraoosane)  for  eaoh  treatment.  For  "heavier" 
eomponents,  whieh  beeame  C-enriehed  during  sliding  exposure,  an  overall  deeease  was 
sometimes  observed  -  however,  eonsiderable  error  (partieularly  with  the  192,000  eyele 
samples)  made  it  diffieult  to  observe  a  model-fitting  trend  (Fig.  5).  Similarly,  some 
lighter  peaks  showing  C-depletion  during  sliding  exposure  inereased  in  eoneentration 
(Fig.  6).  Again,  large  standard  errors  often  made  any  trend  statistieally  insignifieant.  Even 
with  sonieation  and  attempted  full  reeovery,  quantitative  analysis  was  diffieult.  A  PCA 
analysis  did  not  show  any  signifieant  trends  -  due  to  the  relatively  high  replieate 
variability  (data  not  shown).  Beeause  a  C  enriehment  in  some  higher  moleeular  weight 
eomponents  was  observed  with  a  eoneomitant  C  depletion  in  some  lighter  eomponents 
(Figs.  3-4),  and  at  least  in  some  eases  eoneentrations  of  these  eomponents  eo-varied  with 
isotopie  shifts,  we  hypothesize  that  either  terminal  branehed  C-C  or  baekbone  C-C  bonds 
might  have  been  eleaved  during  wear  exposure  -  or  some  eombination  thereof.  A  higher 
degree  of  frietion  stress  and  detailed  eomponent  struetural  analysis  eould  likely  help 
resolve  the  exaet  relationship  between  lubrieant  use  and  shifts  in  eoneentrations  and 
earbon  isotope  ratios.  However,  in  this  study,  analyzing  shifts  in  isotopie  ratios  appeared 
more  robust  in  response  to  frietion  stress  than  monitoring  ohanges  in  eomponent 
eoneentrations. 

3.3.  Comparison  to  "standard"  FTIR  analyses 

FTIR  analysis  has  been  used  to  assess  lubrieant  degradation  and  oontaminant 
aooumulation  [19,  20].  FTIR  analysis  is  relatively  inexpensive,  rapid  and  is  able  to 
disoern  regions  within  a  lubrioants  moleeular  struoture  undergoing  ohange.  FTIR  was 
oonduoted  in  this  study  to  potentially  validate  results  from  isotopie  analysis.  In  order  to 
understand  ohanges  due  to  frietion  stress,  solvent-oorreoted  FTIR  speotra  were  oombined, 
olassed  and  modeled  by  PCA.  The  first  eomponent  (PCI)  aooounted  for  89%  of  the 
speotral  variation  but  only  showed  statistioally-signifioant  variation  (P  <  0.05  with 
frietion  stress  at  384,000  eyeles  (Fig.  7).  To  aooount  for  potential  lubrieant 
oontamination,  two  fuels  were  added,  diesel  and  JP8  (Hero -A  is  an  aviation  lubrieant). 
Fuel  oontaminants  are  often  found  in  used  lubrioants  and  have  been  shown  to  bias  FTIR 
speotra  [20].  Lubrioants  are  generally  higher  in  moleeular  weight  than  fuels,  so  these 
oontaminants  would  not  influenoe  the  GC-C-IRMS  results  (elute  far  earlier  than  lubrieant 
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peaks).  Therefore  eoneomitant  fuel  additions  were  not  analyzed  by  GC-C-IRMS.  The 
fuel-eontaminated  samples  demonstrated  different  variation  with  frietion  stress  -  JP8 
eontaminated  FTIR  speetra  showed  an  initial  inerease  along  PCI  with  a  final  deerease  at 
384,000  cycles.  Diesel-contaminated  lubricant  showed  significant  change  along  PCI  after 
96,000  cycles,  however,  PCI  accounted  for  only  67%  of  the  spectral  variation  (Fig.  7). 

While  the  FTIR  data  did  not  allow  a  robust  model  linking  spectral  features  to  friction 
stress,  the  PCI  loadings  did  show  spectral  changes  occurred  in  the  CH3  and  CH2  regions 
as  well  as  C=0,  other  carbon-oxygen  bonds  and  0-H  bonds.  The  0-H  bond  loading 
appeared  to  vary  significantly  between  neat  lubricant  and  added  fuel  contaminants  (Fig. 

8)  -  showing  no  friction-stress  changes  in  neat  lubricant,  but  perhaps  spurious  changes 
with  fuel-contaminated  lubricant  (Fig.  8;  '-'3200-3400  cm'^).  Aside  from  the  0-H  bond 
assignment,  spectral  variation  as  modeled  in  PCI,  was  within  similar  bond  assignment 
regions  for  each  friction  stress  treatment.  Most  regions  appeared  to  not  be  impacted  by 
fuel  contaminants.  Overall,  the  PCA  model(s)  from  either  neat  or  fuel-contaminated 
lubricant  did  not  offer  robust  linear  or  exponential  relationships  to  friction  stress  (Fig.  7). 
Within  the  neat  samples,  the  variation  seen  primarily  in  C-C  bonds,  the  FTIR  analysis 
supports  results  seen  in  the  stable  isotope  ratios  suggesting  primarily  C-C  bonds  (and  not 
bonds  involving  hydrogen)  were  primary  regions  of  change  due  to  friction  stress  on  the 
lubricant. 

4,  Conclusions 

This  manuscript  explores  whether  isotopic  shifts  might  be  a  viable  alternative  method 
for  assessing  the  friction  stress  “history”  in  lubricants.  We  observed  statistically- 
significant  shifts  in  collective  isotope  ratios  under  relatively  mild  sliding  conditions.  A 
modified  MANOVA  analysis  showed  each  level  of  friction  stress  significantly  altered  the 
collective  carbon  isotope  ratios  of  the  lubricant  components  {P  <  0.01).  A  PCA  model 
was  formulated  that  demonstrated  a  nonlinear  relation  between  5  C  and  friction  stress. 
This  model  was  quite  good  (r  =  0.96)  and  by  analyzing  the  loadings,  showed  several 
lubricant  components  became  C-enriched  while  others  became  C-depleted  in  response 
to  friction  stress.  Changes  in  concentration  were  not  easily  modeled,  while  FTIR  spectra 
showed  no  statistically-significant  change  (but  did  demonstrate  bonds  most  impacted  by 
friction  stress). 

Isotopic  analysis  may  be  an  independent  means  to  assess  lubricant  changes  relative  to 
friction  stress  and  may  offer  another  alternative  in  an  overall  lubricant  management  plan. 
In  order  for  stable  isotope  analysis  to  be  effective,  several  conditions  must  exist,  1)  the 
lubricant  must  be  amenable  to  GC-C-IRMS  analysis  -  having  multiple  individual 
components  resolvable  by  GC  column  coatings  and  run  conditions,  2)  the  lubricant  must 
be  able  to  be  tested  in  its  unused  state  to  obtain  a  baseline  isotopic  signature,  and  3)  each 
new  batch  of  lubricant  must  be  tested  as  isotopic  ratios  vary  based  on  different  base 
feedstocks.  The  overall  advantages  of  the  proposed  method  are  that  it: 

•  does  not  require  accurate  quantitative  recovery  of  lubricants  or  ad  hoc 
development  of  peak  area  ratios; 

•  is  not  sensitive  to  fuel  contaminants  in  the  lubricant  as  components  are  separated 
by  GC  during  the  analysis. 
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•  may  be  more  sensitive  to  small  ehanges  in  lubrieant  eomponents  as  isotopic 
models  related  to  friction  stress  proved  more  robust  and  statistically  significant 
than  peak  area  shifts  or  FTIR  analysis. 

•  does  not  require  an  a  priori  knowledge  of  lubricant  components  -  only  that  they 
can  be  separated  reliably  by  retention  time. 

•  does  not  require  absolute  isotope  calibration  -  the  statistical  tests  need  only  be 
relative  to  a  surrogate  standard  (which  can  be  set  to  an  isotopic  value  of  0). 

The  overall  disadvantages  to  the  method  are 

•  higher  overall  cost  relative  to  FTIR  or  standard  GC  separations 

•  each  lubricant  batch  needs  to  be  fingerprinted  prior  to  use 

•  initial  modeling  could  be  expensive  and  may  lead  to  less  significant  models  for  a 
given  lubricant  (this  paper  reports  on  a  proof  of  concept). 

Tribometric  tests  appropriate  to  the  lubricant’s  use  (z.e.  high  temperature,  high  or  low 
speed,  differing  boundary  conditions,  etc)  would  need  to  be  conducted  for  a  given 
lubricant  at  relevant  conditions  to  vet  the  method  within  a  lubricant  management 
program.  Results  from  these  experiments  could  then  be  modeled  (PCA,  for  instance)  and 
perhaps  scaled  to  process  control  level  chemometrics  so  that  routine  analyses  could  be 
conducted  at  appropriate  timescales.  Tests  involving  lubricant  component  isotope  ratios 
would  be  independent  of  other  protocols  currently  employed  to  assess  lubricant  friction 
stress.  Isotope  evaluation  could  be  coupled  to  lubricant  efficacy  tests  and  may  thus  serve 
as  an  indicator  of  lubricant  senescence.  Additional  experiments  would  need  to  be 
conducted,  probably  using  more  aggressive  friction  stress.  In  addition,  a  lubricant 
management  plan  would  need  to  annotate  initial  sampling  of  all  “new”  lubricant  batches 
for  wear  comparisons.  At  present,  this  proof  of  concept  study  demonstrates  using  isotopic 
ratio  shifts  as  an  independent  method  for  assessing  friction  stress  in  lubricants. 
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Table  1.  Tribometer  run  characteristics 


Cycles 
(X  1000) 

Run  Time  {hr) 

Linear 

Distance  {km) 

0 

0 

0 

24 

4 

2.26 

48 

8 

4.53 

96 

16 

9.05 

192 

32 

18.1 

384 

64 

36.2 
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Table  2.  Probabilities  that  5D  values  represent  the  same  means  among  wear  cycle 
treatments. 


Sample  1 

Sample  2 

P 

C.L 

0  cycles 

24K  cycles 

0.031 

0.0027 

0  cycles 

48K  cycles 

0.85 

0.022 

0  cycles 

96K  cycles 

0.52 

0.027 

0  cycles 

192K  cycles 

0.47 

0.020 

0  cycles 

384K  cycles 

0.14 

0.0037 

24K  cycles 

48K  cycles 

0.61 

0.042 

24K  cycles 

96K  cycles 

0.11 

0.0037 

24K  cycles 

192K  cycles 

0.20 

0.0054 

24K  cycles 

384K  cycles 

0.22 

0.011 

48K  cycles 

96K  cycles 

0.49 

0.011 

48K  cycles 

192K  cycles 

0.32 

0.010 

48K  cycles 

384K  cycles 

0.33 

0.0087 

96K  cycles 

192K  cycles 

0.71 

0.030 

96K  cycles 

384K  cycles 

0.97 

0.11 

192K  cycles 

384K  cycles 

0.93 

0.19 
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Table  3.  Probabilities  that  5  C  values  represent  the  same  means  among  wear  eycle 
treatments. 


Sample  1 

Sample  2 

P 

C.L 

0  eyeles 

24K  eyeles 

0.0005 

0.00001 

0  eyeles 

48K  eyeles 

0.0023 

0.00004 

0  eyeles 

96K  eyeles 

0.0012 

0.00002 

0  eyeles 

192K  eyeles 

0.0076 

0.00010 

0  eyeles 

384K  eyeles 

0.0405 

0.00124 

24K  eyeles 

48K  eyeles 

0.0310 

0.00384 

24K  eyeles 

96K  eyeles 

0.0011 

0.00001 

24K  eyeles 

192K  eyeles 

0.0007 

0.00001 

24K  eyeles 

384K  eyeles 

0.0014 

0.00003 

48K  eyeles 

96K  eyeles 

0.0036 

0.00003 

48K  eyeles 

192K  eyeles 

0.0013 

0.00001 

48K  eyeles 

384K  eyeles 

0.0062 

0.00007 

96K  eyeles 

192K  eyeles 

0.0020 

0.00009 

96K  eyeles 

384K  eyeles 

0.0128 

0.00032 

192K  eyeles 

384K  eyeles 

0.0054 

0.00021 
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Figure  1.  PCI  scores  from  dl3C  PCA  model  related  to  wear  cycles 

Figure  2.  Bi-plot  for  dl3C  PCA  model.  Regions  (A)  and  (B)  show  heavily 
loaded  regions  along  PCI 

Figure  3.5  C  enrichment  related  to  friction  stress  for  various  components 

Figure  4.5  C  depletion  related  to  friction  stress  for  various  components 

Figure  5.  Change  in  relative  peak  area  related  to  friction  stress  for  various  components 

Figure  6.  Change  in  relative  peak  area  related  to  friction  stress  for  various  components 

Figure  7.  FTIR  PCI  scores  related  to  friction  stress 

Figure  8.  PCI  loadings  for  FTIR  spectral  changes  in  response  to  friction  stress 
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Figure  5. 
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